Generalized Hofstadter Model: Band Topology Transitions and Fractional Quantum Hall States 
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Lattice models subjected to an external effective magnetic field can form topological nontrivial bands char- 
acterized by integer Chern numbers. In this paper, we investigate a generalized Hofstadter model with tunable 
nearest-neighbor and next nearest-neighbor hopping, and show that the band topology is much richer compared 
with the conventional Hofstadter model with only the nearest-neighbor hopping. The phase diagram of band 
Chern numbers are obtained numerically for simple rational flux density and a classification of phases is dis- 
cussed. The fractional quantum Hall states in both \C\ — 1 bands and \C\ > 1 bands are used to reflect the 
band topology in different phases. When our model reduces to a one-dimensional lattice, the ground states are 
crucially different from fractional quantum Hall states. 
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Introduction. — Fractional quantum Hall (FQH) states pro- 
vide the most prominent examples of topologically ordered 
phases of matter [1,2]. They support fractionally charged ex- 
citations that are essential resources for topological quantum 
computation [3]. Although FQH states have only been ob- 
served in continuum solid state systems until now, there are 
many proposals of realization of these exotic states for cold 
atoms in optical lattices [4-7]. Numerical simulations of small 
systems also show that various FQH states can indeed exist in 
lattices subject to a uniform perpendicular effective magnetic 
field [5, 8-12]. 

The starting point to investigate the rich FQH physics in 
the presence of an underlying lattice and an effective field is 
the standard Hofstadter model, which normally considers a 
square lattice with complex nearest-neighbor hopping [13]. 
The single-particle spectrum is the well-known Hofstadter 
butterfly with a fractal structure, possessing q bands when the 
magnetic flux density is a rational number p/q (p and q are 
integers which are prime to each other). Each band carries a 
quantized Hall conductance and the band topology can be la- 
beled by an integer Chern number [14, 15]. Because there is 
only nearest-neighbor hopping in the original model, a ques- 
tion arises naturally: can we produce new band topology by 
introducing longer-range hopping? 

In this paper, we systematically visit this problem by con- 
sidering a generalized Hofstadter model on a square lattice 
with tunable nearest-neighbor and next nearest-neighbor hop- 
ping. We numerically calculate the Chern number of each 
band for simple rational <f) = p/q, and find various phases 
of band topology by tuning the strengths of nearest-neighbor 
and next nearest-neighbor hopping for a fixed <\>. The conven- 
tional Hofstadter model corresponds to only one phase in our 
rich phase diagram, and the band topology in other phases are 
strikingly different from it. A classification of phases can be 
given for p = 1 case. In order to identify the band topology 
in a many-body level, we also consider FQH states that can be 
harbored by topological non-trivial bands. By choosing ap- 
propriate hopping parameters, we indeed find different FQH 
states for different band topology. These states are expected to 



be similar to recently discovered fractional Chern insulators 
(FCIs) [16-31] because of the adiabatic continuity between 
Hofstadter and Chern insulator states [32-34]. Finally, we 
study a one-dimensional lattice model for which our general- 
ized Hofstadter model can be regarded as a two-dimensional 
ancestor. We find ground states with crucial difference to FQH 
states at fractional filling factors [35, 36]. 

Model and band topology. — Here we consider a two- 
dimensional (2D) generalized Hofstadter model on a square 
lattice with tunable nearest-neighbor (NN) and next-nearest- 
neighbor (NNN) hopping, 

#2D = - E [< m C„ +1 , m +A°V 2 ^"+5)4 m C„ +1 , m+1 
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where (n, to) is the site coordinate (see Ref. [11] for another 
generalization with arbitrarily long-range hopping exponen- 
tially decayed with the hopping distance). We set t = 1 as the 
energy unit. Because of the complex hopping terms, the 2D 
lattice is pierced by a uniform perpendicular effective mag- 
netic field which breaks the time-reversal symmetry, enabling 
us to label the topological property of an isolated band (with- 
out level crossings with other bands) with eigenstate |k) by 
its Chern number, defined as C — J k BZ f xy (k)d 2 k, where 
f xy (k) = d x a y - d y a x and a,j(k) = -i{k\dj\k) , j = x,y. 

For the conventional Hofstadter model (A od = 0), the 
Chern number of each band can be described by a simple 
picture [15]. In the case of p = 1, the middle band has 
C = — (q — 1) and other bands have C — \ for odd q, while 
the middle two bands have total C = —(q — 2) and other bands 
still have C = 1 for even q. In the case of p > 1, the bands are 
grouped into Q + 1 clusters by writing p/q = 1/(Q +p' /q'), 
and the statement above for p = 1 is still valid for the to- 
tal Chern number in each cluster. However, in the presence 
of the tunable A d and A od , we do not have such a simple pic- 
ture of band topology. By numerically calculating the Chern 
number of each band [37], we find that there are some critical 



values of A od and A d at which the band touching occurs, and 
the Chern numbers of some bands may change after cross- 
ing these critical values, indicating a phase transition of band 
topology. 

In Fig. 1(a), we take <f>= 1/3 as an example to demonstrate 
the rich phase diagram on the A od -A d plane. The whole dia- 
gram is divided into several distinct regions that are symmetric 
with respect to (A od , A d ) = (0, 0). In each region, we can label 
the band topology by Chern numbers of three bands from the 
bottom (with the lowest energy) to the top (with the highest 
energy), i.e. (6i, 6 2 , 6 3 ) satisfying Yn=i Q = 0. The con- 
ventional Hofstadter model corresponds to the axis of A d (ex- 
cept the origin) in phase I, where (C%, 62, 63) = (1, —2, 1) 
is consistent with the conclusion of Ref. [15]. However, the 
band topology in other three phases is strikingly different from 
that in phase I. Chern numbers change due to the band touch- 
ing and a phase transition occurs on the boundary between 
two neighbor regions. The sign effect of A d (A od ) can be seen 
clearly: changing the sign of A d (A od ) induces a flip of Chern 
numbers from (61,62,63) to (6*3,(72, Ci). We can clas- 
sify all phases into two classes according to their Chern num- 
bers. In the first class (phase I, II and III), (Ci, 62,63) = 
V{1, —2, 1) (P means permutation), and in the second class 
(phase IV), ((71,(72,(73) = (-2,4,-2). 

For other values of <f> with larger q, similar but more compli- 
cated phase diagrams can be observed. In Fig. 1(b), we show 
the phase diagram for (\> = 1/4. There are six phases differ- 
entiated by Chern numbers of four bands from the bottom to 
the top, i.e. (61 ,6*2,(73, 64), and the conventional Hofstadter 
model is located on the boundary between phase I and phase 
II. Similar to the <p = 1/3 case, the diagram is symmetric 
with respect to (A od , A d ) = (0, 0) and the Chern numbers are 
flipped by the sign change of A d or A od . All phases can also 
be classified into two classes. In the first class (phase I, II, V 
and VI), (61, 6 2 , 63, 64) = V(l, 1, -3, 1), and in the second 
class (phase III and IV), (61, 6 2 , 63, 6 4 ) = 7^(1, -3,5,-3 ) 
[the line under (—3, 5, —3) means that they are grouped into a 
cluster and should be moved as a whole in the permutation]. 

The symmetry and the classification of phases in the 
phase diagram that we observe for tp — 1/3 and cf> = 1/4 
are inherited by <p = 1/5 [38]. Therefore, we expect 
that there are 2q — 2 phases which can be classified 
into two classes (i) and (ii) for <fi = 1/q. We have 

(61,62, ..., 6|-g/ 2 JiC'[g/2jH-l5 C*rg/2j+2i ■■■,C q -i,C q ) 

V(l, 1, ..., 1, 1 — q, 1, ..., 1, 1) in class (i) with q phases, while 
(61,62, ..., C\ q /2\,C\ q /2\+\, 6p g /2j + 2, ...,6 g _i,6 9 ) 
V(l, 1, •■-, 1— q, 1 + <Z, 1— q, • ■•, 1, 1) in class (ii) with q — 2 
phases, where \x\ means the integer part of x. For </> = p/q 
with p > 1, the situation is much more complicated. How- 
ever, most of the (A od , A d ) plane is occupied by phases in 
which the band Chern numbers (61, ..., C q ) can be obtained 
by the permutation of that for the conventional Hofstadter 
model [38]. 

Fractional quantum Hall states. — In order to further char- 
acterize the band topology of different phases discovered in 
the last section, we consider interacting particles partially 




Figure 1 . (Color online) The phase diagram in the A od -A d plane for 
(a) cf> — 1/3 and (b) cj> — 1/4. (a) There are four phases with 
different band Chern numbers: (I) (d, C 2 , 63) = (1,-2,1); (II) 
(Ci, C 2 ,6 3 ) = (1,1,-2); (III) (61,62,63) = (-2,1,1); (IV) 
(61,62,63) = (-2,4,-2). (b) There are six phases with dif- 
ferent band Chern numbers: (I) (61,62,63,64) = (1,1,-3,1); 
(II) (61,62,63,64) = (1,-3,1,1); (III) (61,62,63,64) = 
(1,-3,5,-3); (IV) (61,62,63,64) = (-3,5,-3,1); (V) 
(61,62,63,64) = (1,1,1,-3); (VI) (61,62,63,64) = 
(—3, 1, 1, 1). On the boundary between two regions, band touching 



filled in one band. If the Chern number of this occupied band 
changes due to a phase transition, the FQH states that this 
band favours to harbor should also change. Therefore, we can 
utilize FQH states to reflect the band topology. In the fol- 
lowing, we choose <\> = 1/3 to study the FQH states in var- 
ious phases. These states are expected to be similar to FCIs 
in topological flat bands [16-31], although the net magnetic 
field is nonzero in our model (a gauge transformation can be 
used to obtain a zero net magnetic field [32]). 

We first consider N e fermions partially filled in the mid- 
dle band with NN and NNN repulsive interaction H m = 

V i H{i,j) ninj+V 2 E<(ij» ™* n j' with (*>i) and «*» J» rep- 
resenting NN and NNN sites respectively. In order to focus 

on the topological property of the band, we take the flat band 
limit and diagonalize H mt projected to the flattened occupied 
band for a finite system on the torus with iVi x N2 unit cells. 
We choose three sites in the x direction as a unit cell, so the 
actual lattice size is 3iVi x JVjj. The filling factor v is defined 
as N e /(NiN2). Each energy level can be labeled by the 2D 
total momentum (Ki, Kq). 

In phase II, the Chern number of the middle band is 62 = 1, 
so v = 1/3 fermionic FQH state may be stabilized for appro- 
priate hopping parameters. By choosing (A d , A od ) = (1,0.5), 
we indeed find three quasidegenerate ground states at the bot- 
tom of the energy spectrum separated by an energy gap from 
high excited levels [Fig. 2(a)]. The spectral flow under twisted 
boundary conditions also confirms that the ground states are 
topologically nontrivial [Fig. 2(b)]. When the boundary phase 
$ changes from to 3 x 27r, the three ground states evolve 
into each other, being always separated from excited states by 
a gap, and finally return to the initial configuration. In order to 
discard competing possibilities, such as charge density waves, 
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Figure 2. (Color online) Numerical results for fermionic v = 1/3 
[(a-c)] and v = 1/5 [(d-f)] FQH states with Vi = V 2 = 0.5. 
(A d ,A od ) = (1,0.5) for (a-c) and (A d , A od ) = (0.75,0) for (d-f). 
(a) The low-energy spectra for iV e = 6, 8, 10. (b) The x-direction 
spectral flow for iV e = 8. (c) The PES for N e = 9, N A = 4, and 
JVi x N% = 3x9. The number of states below the gap is 5508. 
(d) The low-energy spectra for N e = 5, 6, 7. (e) The y-direction 
spectral flow for N e =6. (f) The PES for N e = 7, N A = 3 and 
Ni x JV2 = 5x7. The number of states below the gap is 2695. 
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Figure 3. (Color online) Numerical results for bosonic v = 1/3 
[(a-c)] and v = 1/2 [(d-f)] FQH states with U = 2V X = 1.0. 
(A d ,A od ) = (-0.75,5) for (a-c) and (A d ,A od ) = (1,0.5) for (d- 
f). (a) The low-energy spectra for Nb = 4, 6, 8. (b) The a;-direction 
spectral flow for N b = 5. (c) The PES for N b = 8, Ni x N 2 = 3 x 8 
and Na = 4. The number of states below the gap is 2730. (d) The 
low-energy spectra for JVj, = 4, 6, 8. (e) The y-direction spectral 
flow for N b = 6. (f) The PES for N b = 9, iVi x iV 2 = 3 x 6 and 
Na = 4. The number of states below the gap is 1287. 



we investigate the particle-cut entanglement spectrum (PES) 
[20, 30, 39, 40], which can probe the excitation structure of 
the ground states [Fig. 2(c)]. One can see a clear entangle- 
ment gap in the spectrum, and the number of low-lying levels 
below the gap exactly matches the quasihole counting in the 
corresponding Abelian FQH states. Combining all these evi- 
dences together, we are convinced that the v — 1/3 fermionic 
FQH state exists in the middle band in phase II. 

After the system evolves from phase II to phase I, the Chern 
number of the middle band changes from C-2 = 1 to C<i = —2. 
Therefore, it is expected to observe the v = 1/5 instead of the 
v = 1/3 fermionic FQH state for appropriate hopping param- 
eters [30]. Our numerical results for (A d , A od ) = (0.75,0), in- 
cluding the energy spectrum, spectra flow, and the PES, pro- 
vide convincing evidences of the v = 1/5 fermionic FQH 
state [Fig. 2(d-f)] in C = -2 band. Moreover, a v = 1/3 
state in the middle band like that appearing in phase II is not 
found in phase I. In this way, we characterize the different 
band topology between phase I and phase II in a many-body 
level. 

Now we turn our attention to N b bosons partially filled 
in the lowest band with on-site and NN repulsive interaction 

#int = Uj2t n ii n i ~ 1) + V52(ij) n i n j- In P haSe IV > me 

Chern number of the lowest band is C\ = —2, which may 
harbor the unusual bosonic FQH state at v = 1/3 with an odd 



denominator [29, 30]. Our numerical results indeed support 
the existence of this state [Fig. 3(a-c)]. Moreover, in phase II 
where C\ changes to 1 from —2, we find that bosons form a 
FQH state at v = 1/2 instead of v = 1/3 [Fig. 3(d-f)]. There- 
fore, the different band topology between phase II and phase 
IV is also confirmed in a many-body picture. 

Interacting one -dimensional model. — Our 2D generalized 
Hofstadter model can be regarded as an ancestor of a one- 
dimensional (ID) model with the NN hopping and the on-site 
potential governed by a modulation of frequency <j>, 

H 1D = - Y^ { [1 + 2A od cos(2n(j)n + 5 + n<f>)] c{c n+1 + h.c 

n 

+2A d cos(2tt<M + <*)4cn I , (2) 

where 6 is a phase factor that can be randomly chosen. The 
topological properties in the case of A od = were stud- 
ied in Ref. [41]. After replacing c„ with c n .s in Eq. (2), 
we can return to our generalized 2D Hofstadter model by 

#2D = h Jo' H 1D (6)dS. 

Motivated by recent discoveries of gapped phases of in- 
teracting particles in ID flat bands at fractional fillings [36, 
42, 43], we consider the lattice model in Eq. (2) under peri- 
odic boundary condition partially filled by N e fermions with 
dipole-dipole interaction H mt = J2i^j ul^-U f° r 4> = 1/3- 
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Figure 4. (Color online) Numerical results for the ID model in 
Eq. (2) at v = 1/3 [(a-b)] and v = 1/5 [(c-d)], with S = tt/2, 
(j> = 1/3, and (A d , A od ) = (0.75, 0). The interaction is projected to 
the middle band, (a) The low-energy spectrum at v — 1/3. (b) The 
PES for N e = 7, iVceii = 21 and A^ = 3. The number of states 
below the gap is 105. (c) The low-energy spectrum at v — 1/5. (d) 
The PES for N e = 5, JV c di = 25 and N A = 2. The number of states 
below the gap is 50. 



The filling factor v is defined as N e /N ce n with N ce \\ the num- 
ber of unit cells, each of which contains three sites. In the 
flat band limit, we find that there are always m-fold ground- 
state degeneracy at v = 1/m, independent of the values of 
(A d , A od ) (Fig. 4). Moreover, the number of states below the 
entanglement gap in the PES is much smaller than the FQH 



'N A 



mN e - (m - 1)N A - 1 
N A -l 



but 



counting W^ = 

matches the charge density wave (CDW) counting A/" cr j^, = 

( N \ 
TO I AJ - e I [35, 36]. Therefore, we conclude that the many- 
body ground states cannot differentiate the band topology in 
the ID model in Eq. (2). 

Conclusions. — In summary, we investigate a generalized 
Hofstadter model by considering tunable nearest-neighbor 
and next-nearest-neighbor hopping on a square lattice sub- 
jected to a uniform perpendicular effective magnetic field. 
Rich phase diagrams of the band topology exist on the A od - 
A d plane for rational flux density <fi — p/q and a classification 
of phases is discussed foip = 1, The many-body FQH states 
can also be used to differentiate the band topology in various 
phases. However, the situation is completely different in ID, 
where the many-body ground states have the CDW property 
and cannot reflect the band topology. It is likely that our work 
would be helpful for the discovery of topological phases of 
matter in cold-atom experiments. 

Acknowledgements. — Z. L. acknowledges Emil Bergholtz 
and Andreas Lauchli for related collaborations and thanks 
Emil Bergholtz for useful discussions. This work is supported 
by "973" program (2010CB922904), NSFC and grants from 
CAS. Z. L. is supported by the China Postdoctoral Science 
Foundation Grant No. 2012M520149. 



* zliu@csrc.ac.cn 
' hfan@iphy.ac.cn 
[1] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
[2] G. Moore, and N. Read, Nucl. Phys. B 360, 362 (1991). 
[3] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das 

Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
[4] D. Jaksch and P. Zoller, New J. Phys. 5, 56 (2003). 
[5] A. S. S0rensen, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 

94, 086803 (2005). 
[6] F. Gerbier and J. Dalibard, New J. Phys. 12, 033007 (2010). 
[7] N. R. Cooper, Phys. Rev. Lett. 106, 175301 (2011). 
[8] R. N. Palmer and D. Jaksch, Phys. Rev. Lett. 96, 180407 (2006). 
[9] R. N. Palmer, A. Klein, and D. Jaksch, Phys. Rev. A 78, 

013609 (2008). 
[10] G. Moller and N. R. Cooper, Phys. Rev. Lett. 103, 105303 

(2009). 
[11] E. Kapit and E. Mueller, Phys. Rev. Lett. 105, 215303 (2010). 
[12] L. Hormozi, G Moller, and S. H. Simon, Phys. Rev. Lett. 108, 

256809 (2012). 
[13] D. R. Hofstadter, Phys. rev. B 14, 2239 (1976). 
[14] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den 

Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[15] M. Kohmoto, Phys. Rev. B 39, 11943 (1989). 
[16] E. Tang, J.-W. Mei, and X.-G. Wen, Phys. Rev. Lett. 106, 

236802(2011). 
[17] K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, 

Phys. Rev. Lett. 106, 236803 (201 1). 
[18] T Neupert, L. Santos, C. Chamon, and C. Mudry, 

Phys. Rev. Lett. 106, 236804 (201 1). 
[19] D. N. Sheng, Z. Gu, K. Sun, and L. Sheng, Nat. Commun. 2, 

389(2011). 
[20] N. Regnault and B. A. Bernevig, Phys. Rev. X 1, 021014 

(2011); B. A. Bernevig and N. Regnault, Phys. Rev. B 85, 

075128 (2012); Y.-L. Wu, B. A. Bernevig, N. Regnault, Phys. 

Rev. B 85, 0751 16 (2012); T Liu, C. Repellin, B. A. Bernevig, 

and N. Regnault, arXiv: 1206.2626; Y.-L. Wu, B. A. Bernevig, 

N. Regnault, Phys. Rev. B 86, 085129 (2012). 
[21] J. W. F. Venderbos, S. Kourtis, J. van den Brink, and 

M. Daghofer, Phys. Rev. Lett. 108, 126405 (2012). 
[22] N. Y Yao, A. V. Gorshkov, C. R. Laumann, A. M. Lauchli, 

Y Ye, and M. D. Lukin, arXiv: 1212.4839. 
[23] N. R. Cooper and R. Moessner, Phys. Rev. Lett. 109, 215302 

(2012). 
[24] A. Lauchli, Z. Liu, E. J. Bergholtz, and R. Moessner, 

arXiv: 1207.6094. 
[25] N. R. Cooper and J. Dalibard, arXiv: 1212.3552. 
[26] Y.-F. Wang, Z.-C. Gu, C.-D. Gong, and D. N. Sheng, Phys. Rev. 

Lett. 107, 146803 (2011); Y.-F. Wang, H. Yao, Z.-C. Gu, C- 

D. Gong, D. N. Sheng, Phys. Rev. Lett. 108, 126805 (2012). 
[27] X.-L. Qi, Phys. Rev. Lett. 107, 126803 (2011). 
[28] M. Trescher and E. J. Bergholtz, Phys. Rev. B 86, 241111(R) 

(2012). 
[29] Y.-F. Wang, H. Yao, C.-D. Gong, and D. N. Sheng, Phys. Rev. 

B86, 201 101(R) (2012). 
[30] Z. Liu, E. J. Bergholtz, H. Fan, and A. M. Lauchli, Phys. Rev. 

Lett. 109, 186805 (2012). 
[31] A. Sterdyniak, C. Repellin, B. A. Bernevig, and N. Regnault, 

arXiv: 1207.6385; Y.-L. Wu, N. Regnault, and B. A. Bernevig, 

Phys. Rev. Lett. 110, 106802 (2013). 
[32] Y-H. Wu, J. K. Jain, and K Sun, Phys. Rev. B 86, 165129 

(2012). 



[33] Z. Liu and E. J. Bergholtz, Phys. Rev. B 87, 035306 (2013). 

[34] T. Scaffidi, and G. Moller, Phys. Rev. Lett. 109, 246805 (2012). 

[35] B. A. Bernevig and N. Regnault, arXiv: 1204.5682. 

[36] J. C. Budich and E. Ardonne, arXiv: 1304.4366. 

[37] T. Fukui, Y. Hatsugai, and H. Suzuki, J. Phys. Soc. Jpn. 74, 

1674 (2005). 
[38] See the Supplementary Material for the phase diagrams for <f> = 

1/5 and <j> = 2/5. 
[39] A. Sterdyniak, N. Regnault, and B.A. Bernevig, Phys. Rev. Lett. 



106, 100405(2011). 
[40] A. Chandran, M. Hermanns, N. Regnault, and B.A. Bernevig, 

Phys. Rev. B 84, 205136 (2011). 
[41] L.-J. Lang, X.-M. Cai, and S. Chen, Phys. Rev. Lett. 108, 

220401 (2012). 
[42] H.-M. Guo, S.-Q. Shen, S.-P Feng, Phys. Rev. B 86, 085124 

(2012); H.-M. Guo, Phys. Rev. A 86, 055604 (2012). 
[43] Z.-H. Xu, L.-H. Li, and S. Chen, arXiv: 1210.7696. 



SUPPLEMENTARY MATERIAL FOR "GENERALIZED HOFSTADTER MODEL: BAND TOPOLOGY TRANSITIONS AND 

FRACTIONAL QUANTUM HALL STATES" 

In the main text, we have shown the rich phase diagrams for = 1/3 and <fi = 1/4. The number of distinct phases, where the 
nontrivial band topology is labeled by the Chern numbers, increases as the denominator of the rational flux density (</> = p/q) 
becomes larger. Here, we will show the results for <p = 1/5 and <fi = 2/5, where more phases exist than <j> = 1/3 and <fi = 1/4 
cases. Since changing the sign of A d or A od only flips the band Chern numbers, we only illustrate the phase diagram in the first 
quadrant of A d -A od plane: (A d > 0, A od > 0). 



o 




c 




Figure 5. (Color online) (a) The phase diagram for <j> = 1/5, with only the first quadrant (A d > 0, A od > 0) being presented. Phase I- VI can be 
seen. Phase VII and VIII are located in the second quadrant, (b) The phase diagram for <f> = 2/5, with only the first quadrant (A d > 0, A od > 0) 
being presented. Phase I- VI can be seen. Phase VII and VIII are located in the second quadrant. Phase IX is located beyond A od = 5. 



The phase diagram for <fi — 1/5 is shown in Fig. 5(a). Six of the eight different phases represented by different colors are il- 
lustrated in the diagram, and other two absent phases are located in the second quadrant of the A d -A od plane: (A d < 0, A od > 0). 
The band Chern numbers are numerically calculated in each phase, and the result is exhibited in Table I, The conventional 
Hofstadter model is located in phase I, We can see that for phase I and IV, the band Chern number (C\, C2, C3, C4, C5) is 
symmetric, i.e. C\ = C5 and C2 = C4. That is because, by taking the diagonal limit (A d 7^ 0, A od = 0) or the off-diagonal 
limit (A d = 0, A od 7^ 0), the Hamiltonian preserves the particle-hole symmetry, which makes the band Chern numbers sym- 
metric. Furthermore, according to the data in Table I, we can see that the eight phases can be divided into two classes: (i) 
(Ci, C* 2 , C 3 , C 4 , C 5 ) = V(l, 1, -4, 1, 1) and (ii) (Ci, C 2 ,C 3 , C 4 , C 5 ) = V(l, -4,6,-4 , 1), which is similar to the = 1/3 
and 4> = 1/4 cases. The classification of phases for a general <j> = l/q is given in the main text. 

However, the situation is different forp^l. Here we take <f> = 2/5 as an example to demonstrate its different behavior from 
6 = 1/5. The phase diagram is shown in Fig. 5(b), and the Chern number of each band is listed in Table II. The conventional 
Hofstadter model (A od = 0) is located in phase I, where the band Chern numbers (C\, C2, C3, C4, C5) = (—2, 3, —2, 3, —2). 



This is consistent with the conclusion that the five bands should be grouped into three clusters with the total Chern num- 
ber 1, —2 and 1, respectively. We write {C\, C<i, C3, C4, C5) = (—2,3,-2,3,-2) to show such a clustering picture. The 
nine phases in Table II can be divided into three classes: (i) (Ci, C2, C3, C4, C5) = V{— 2, 3,— 2, 3, — 2), i.e. phase I-V 
and VII; (ii) (C*i, C 2 , C 3 , C 4 , C* 5 ) = (-2, 3, 3, -7, 3) and its inversion, i.e. phase VI and VmT(iii) (C^C 2 , C 3 , C 4 , C 5 ) = 
(—2, —2, 8, —2, —2), i.e. phase IX. One can notice that most of the phase space is occupied by phases belonging to class (i). 

Now, we can have a more complete look at how the conventional Hofstadter model is generalized to our model. The anisotropic 
hopping and NNN hopping lead to the change of band structure and band topology and more topological nontrivial phases 
emerge. Most of these new phases can be regarded as daughters of the conventional Hofstadter phase in the sense that their band 
Chern numbers are just permutations of that in the father model. 

Table I. The band Chern numbers in all phases for cf> — 1/5. In each column, the Chern numbers are listed from band 5 on the top (with the 
highest energy) to band 1 at the bottom (with the lowest energy). 
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Table II. The band Chern numbers in all phases for <f> — 2/5. In each column, the Chern numbers are listed from band 5 on the top (with the 
highest energy) to band 1 at the bottom (with the lowest energy). 
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